Computational methods in the analysis of SARS-CoV-2 in mammals: A systematic review of the literature.

SARS-CoV-2 is an enveloped RNA virus that causes severe respiratory illness in humans and animals. It infects cells by binding the Spike protein to the host's angiotensin-converting enzyme 2 (ACE2). The bat is considered the natural host of the virus, and zoonotic transmission is a significant risk and can happen when humans come into close contact with infected animals. Therefore, understanding the interconnection between human, animal, and environmental health is important to prevent and control future coronavirus outbreaks. This work aimed to systematically review the literature to identify characteristics that make mammals suitable virus transmitters and raise the main computational methods used to evaluate SARS-CoV-2 in mammals. Based on this review, it was possible to identify the main factors related to transmissions mentioned in the literature, such as the expression of ACE2 and proximity to humans, in addition to identifying the computational methods used for its study, such as Machine Learning, Molecular Modeling, Computational Simulation, between others. The findings of the work contribute to the prevention and control of future outbreaks, provide information on transmission factors, and highlight the importance of advanced computational methods in the study of infectious diseases that allow a deeper understanding of transmission patterns and can help in the development of more effective control and intervention strategies.

Evolution of human H3N2 influenza virus receptor specificity has substantially expanded the receptor-binding domain site.

Hemagglutinins (HAs) from human influenza viruses descend from avian progenitors that bind α2-3-linked sialosides and must adapt to glycans with α2-6-linked sialic acids on human airway cells to transmit within the human population. Since their introduction during the 1968 pandemic, H3N2 viruses have evolved over the past five decades to preferentially recognize human α2-6-sialoside receptors that are elongated through addition of poly-LacNAc. We show that more recent H3N2 viruses now make increasingly complex interactions with elongated receptors while continuously selecting for strains maintaining this phenotype. This change in receptor engagement is accompanied by an extension of the traditional receptor-binding site to include residues in key antigenic sites on the surface of HA trimers. These results help explain the propensity for selection of antigenic variants, leading to vaccine mismatching, when H3N2 viruses are propagated in chicken eggs or cells that do not contain such receptors.

Structural insights into the interaction between adenovirus C5 hexon and human lactoferrin.

Adenovirus (AdV) infection of the respiratory epithelium is common but poorly understood. Human AdV species C types, such as HAdV-C5, utilize the Coxsackie-adenovirus receptor (CAR) for attachment and subsequently integrins for entry. CAR and integrins are however located deep within the tight junctions in the mucosa where they would not be easily accessible. Recently, a model for CAR-independent AdV entry was proposed. In this model, human lactoferrin (hLF), an innate immune protein, aids the viral uptake into epithelial cells by mediating interactions between the major capsid protein, hexon, and yet unknown host cellular receptor(s). However, a detailed understanding of the molecular interactions driving this mechanism is lacking. Here, we present a new cryo-EM structure of HAdV-5C hexon at high resolution alongside a hybrid structure of HAdV-5C hexon complexed with human lactoferrin (hLF). These structures reveal the molecular determinants of the interaction between hLF and HAdV-C5 hexon. hLF engages hexon primarily via its N-terminal lactoferricin (Lfcin) region, interacting with hexon's hypervariable region 1 (HVR-1). Mutational analyses pinpoint critical Lfcin contacts and also identify additional regions within hLF that critically contribute to hexon binding. Our study sheds more light on the intricate mechanism by which HAdV-C5 utilizes soluble hLF/Lfcin for cellular entry. These findings hold promise for advancing gene therapy applications and inform vaccine development. IMPORTANCE: Our study delves into the structural aspects of adenovirus (AdV) infections, specifically HAdV-C5 in the respiratory epithelium. It uncovers the molecular details of a novel pathway where human lactoferrin (hLF) interacts with the major capsid protein, hexon, facilitating viral entry, and bypassing traditional receptors such as CAR and integrins. The study's cryo-EM structures reveal how hLF engages hexon, primarily through its N-terminal lactoferricin (Lfcin) region and hexon's hypervariable region 1 (HVR-1). Mutational analyses identify critical Lfcin contacts and other regions within hLF vital for hexon binding. This structural insight sheds light on HAdV-C5's mechanism of utilizing soluble hLF/Lfcin for cellular entry, holding promise for gene therapy and vaccine development advancements in adenovirus research.

Variation in the affinity of three representative avian adenoviruses for the cellular coxsackievirus and adenovirus receptor.

According to previous studies, three representative avian adenoviral strains utilize coxsackievirus-adenovirus receptor (CAR) as a receptor and seem to exhibit diverse binding affinities and modes. Thus, further revealing the exact molecular mechanism underlying the interaction between different FAdVs and the attachment receptor CAR is necessary. In this study, we successfully solved the crystal structure of the FAdV-4 fiber1 knob at 1.6 Šresolution. The interaction between the fibre knob and different domains of CAR was verified by confocal microscopy, coimmunoprecipitation and surface plasmon resonance (SPR) analysis. The fibre knobs of the three representative fowl adenoviruses specifically recognized CAR domain 1 (D1), but the recognition of CAR domain 2 (D2) by chicken embryo lethal orphan (CELO) strains was weak. These results provide insights into the differences in adenovirus‒host cell interactions and have important implications for the exploration of viral invasion mechanisms.

Predicting virus Fitness: Towards a structure-based computational model.

Predicting the impact of new emerging virus mutations is of major interest in surveillance and for understanding the evolutionary forces of the pathogens. The SARS-CoV-2 surface spike-protein (S-protein) binds to human ACE2 receptors as a critical step in host cell infection. At the same time, S-protein binding to human antibodies neutralizes the virus and prevents interaction with ACE2. Here we combine these two binding properties in a simple virus fitness model, using structure-based computation of all possible mutation effects averaged over 10 ACE2 complexes and 10 antibody complexes of the S-protein (∼380,000 computed mutations), and validated the approach against diverse experimental binding/escape data of ACE2 and antibodies. The ACE2-antibody selectivity change caused by mutation (i.e., the differential change in binding to ACE2 vs. immunity-inducing antibodies) is proposed to be a key metric of fitness model, enabling systematic error cancelation when evaluated. In this model, new mutations become fixated if they increase the selective binding to ACE2 relative to circulating antibodies, assuming that both are present in the host in a competitive binding situation. We use this model to categorize viral mutations that may best reach ACE2 before being captured by antibodies. Our model may aid the understanding of variant-specific vaccines and molecular mechanisms of viral evolution in the context of a human host.

Biophysical evolution of the receptor-binding domains of SARS-CoVs.

With hundreds of coronaviruses (CoVs) identified in bats that can infect humans, it is essential to understand how CoVs that affected the human population have evolved. Seven known CoVs have infected humans, of which three CoVs caused severe disease with high mortalities: severe acute respiratory syndrome (SARS)-CoV emerged in 2002, Middle East respiratory syndrome-CoV in 2012, and SARS-CoV-2 in 2019. SARS-CoV and SARS-CoV-2 belong to the same family, follow the same receptor pathway, and use their receptor-binding domain (RBD) of spike protein to bind to the angiotensin-converting enzyme 2 (ACE2) receptor on the human epithelial cell surface. The sequence of the two RBDs is divergent, especially in the receptor-binding motif that directly interacts with ACE2. We probed the biophysical differences between the two RBDs in terms of their structure, stability, aggregation, and function. Since RBD is being explored as an antigen in protein subunit vaccines against CoVs, determining these biophysical properties will also aid in developing stable protein subunit vaccines. Our results show that, despite RBDs having a similar three-dimensional structure, they differ in their thermodynamic stability. RBD of SARS-CoV-2 is significantly less stable than that of SARS-CoV. Correspondingly, SARS-CoV-2 RBD shows a higher aggregation propensity. Regarding binding to ACE2, less stable SARS-CoV-2 RBD binds with a higher affinity than more stable SARS-CoV RBD. In addition, SARS-CoV-2 RBD is more homogenous in terms of its binding stoichiometry toward ACE2 compared to SARS-CoV RBD. These results indicate that SARS-CoV-2 RBD differs from SARS-CoV RBD in terms of its stability, aggregation, and function, possibly originating from the diverse receptor-binding motifs. Higher aggregation propensity and decreased stability of SARS-CoV-2 RBD warrant further optimization of protein subunit vaccines that use RBD as an antigen by inserting stabilizing mutations or formulation screening.

Receptor Density-Dependent Motility of Influenza Virus Particles on Surface Gradients.

Influenza viruses can move across the surface of host cells while interacting with their glycocalyx. This motility may assist in finding or forming locations for cell entry and thereby promote cellular uptake. Because the binding to and cleavage of cell surface receptors forms the driving force for the process, the surface-bound motility of influenza is expected to be dependent on the receptor density. Surface gradients with gradually varying receptor densities are thus a valuable tool to study binding and motility processes of influenza and can function as a mimic for local receptor density variations at the glycocalyx that may steer the directionality of a virus particle in finding the proper site of uptake. We have tracked individual influenza virus particles moving over surfaces with receptor density gradients. We analyzed the extracted virus tracks first at a general level to verify neuraminidase activity and subsequently with increasing detail to quantify the receptor density-dependent behavior on the level of individual virus particles. While a directional bias was not observed, most likely due to limitations of the steepness of the surface gradient, the surface mobility and the probability of sticking were found to be significantly dependent on receptor density. A combination of high surface mobility and high dissociation probability of influenza was observed at low receptor densities, while the opposite occurred at higher receptor densities. These properties result in an effective mechanism for finding high-receptor density patches, which are believed to be a key feature of potential locations for cell entry.

Exogenous Chemicals Impact Virus Receptor Gene Transcription: Insights from Deep Learning.

Despite the fact that coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been disrupting human life and health worldwide since the outbreak in late 2019, the impact of exogenous substance exposure on the viral infection remains unclear. It is well-known that, during viral infection, organism receptors play a significant role in mediating the entry of viruses to enter host cells. A major receptor of SARS-CoV-2 is the angiotensin-converting enzyme 2 (ACE2). This study proposes a deep learning model based on the graph convolutional network (GCN) that enables, for the first time, the prediction of exogenous substances that affect the transcriptional expression of the ACE2 gene. It outperforms other machine learning models, achieving an area under receiver operating characteristic curve (AUROC) of 0.712 and 0.703 on the validation and internal test set, respectively. In addition, quantitative polymerase chain reaction (qPCR) experiments provided additional supporting evidence for indoor air pollutants identified by the GCN model. More broadly, the proposed methodology can be applied to predict the effect of environmental chemicals on the gene transcription of other virus receptors as well. In contrast to typical deep learning models that are of black box nature, we further highlight the interpretability of the proposed GCN model and how it facilitates deeper understanding of gene change at the structural level.

Biomimetic SARS-CoV-2 Spike Protein Nanoparticles.

COVID-19 is an infectious respiratory disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This virus contains a crucial coat protein that engages with target cells via a receptor binding domain (RBD) on its spike protein. To better study the RBD and its therapeutic opportunities, we genetically engineered a simple fusion with a thermo-responsive elastin-like polypeptide (ELP). These fusions express in Escherichia coli at a high yield in the soluble fraction and were easily purified using ELP-mediated phase separation (79 mg/L culture). Interestingly, they assembled peptide-based nanoparticles (Rh = 71.4 nm), which was attributed to oligomerization of RBDs (25.3 kDa) counterbalanced by steric stabilization by a soluble ELP (73.4 kDa). To investigate their biophysical properties, we explored the size, shape, and binding affinity for the human angiotensin-converting enzyme 2 (hACE2) and cellular uptake. Biomimetic nanoparticles such as these may enable future strategies to target the same cells, tissues, and cell-surface receptors as those harnessed by SARS-CoV-2.

SARS-CoV-2 triggers pericyte-mediated cerebral capillary constriction.

The SARS-CoV-2 receptor, ACE2, is found on pericytes, contractile cells enwrapping capillaries that regulate brain, heart and kidney blood flow. ACE2 converts vasoconstricting angiotensin II into vasodilating angiotensin-(1-7). In brain slices from hamster, which has an ACE2 sequence similar to human ACE2, angiotensin II evoked a small pericyte-mediated capillary constriction via AT1 receptors, but evoked a large constriction when the SARS-CoV-2 receptor binding domain (RBD, original Wuhan variant) was present. A mutated non-binding RBD did not potentiate constriction. A similar RBD-potentiated capillary constriction occurred in human cortical slices, and was evoked in hamster brain slices by pseudotyped virions expressing SARS-CoV-2 spike protein. This constriction reflects an RBD-induced decrease in the conversion of angiotensin II to angiotensin-(1-7) mediated by removal of ACE2 from the cell surface membrane and was mimicked by blocking ACE2. The clinically used drug losartan inhibited the RBD-potentiated constriction. Thus, AT1 receptor blockers could be protective in COVID-19 by preventing pericyte-mediated blood flow reductions in the brain, and perhaps the heart and kidney.

Excavating phytochemicals from plants possessing antiviral activities for identifying SARS-CoV hemagglutinin-esterase inhibitors by diligent computational workflow.

Coronaviruses (CoVs) belong to a group of RNA viruses that cause diseases in vertebrates including. Newer and deadlier than SARS CoV-2 are sought to appear in future for which the scientific community must be prepared with the strategies for their control. Spike protein (S-protein) of all the CoVs require angiotensin-converting enzyme2 (ACE2), while CoVs also require hemagglutinin-acetylesterase (HE) glycoprotein receptor to simultaneously interact with O-acetylated sialic acids on host cells, both these interactions enable viral particle to enter host cell leading to its infection. Target inhibition of viral S-protein and HE glycoprotein receptor can lead to a development of therapy against the SARS CoV-2. The proposition is to recognize molecules from the bundle of phytochemicals of medicinal plants known to possess antiviral potentials as a lead that could interact and mask the active site of, HE glycoprotein which would ideally bind to O-acetylated sialic acids on human host cells. Such molecules can be addressed as 'HE glycoprotein blockers'. A library of 110 phytochemicals from Withania somnifera, Asparagus racemosus, Zinziber officinalis, Allium sativum, Curcuma longa and Adhatoda vasica was constructed and was used under present study. In silico analysis was employed with plant-derived phytochemicals. The molecular docking, molecular dynamics simulations over the scale of 1000 ns (1 µs) and ADMET prediction revealed that the Withania somnifera (ashwagandha) and Asparagus racemosus (shatavari) plants possessed various steroidal saponins and alkaloids which could potentially inhibit the COVID-19 virus and even other CoVs targeted HE glycoprotein receptor.Communicated by Ramaswamy H. Sarma.

Mapping of functional SARS-CoV-2 receptors in human lungs establishes differences in variant binding and SLC1A5 as a viral entry modulator of hACE2.

BACKGROUND: The COVID-19 pandemic is an infectious disease caused by SARS-CoV-2. The first step of SARS-CoV-2 infection is the recognition of angiotensin-converting enzyme 2 (ACE2) receptors by the receptor-binding domain (RBD) of the viral Spike (S) glycoprotein. Although the molecular and structural bases of the SARS-CoV-2-RBD/hACE2 interaction have been thoroughly investigated in vitro, the relationship between hACE2 expression and in vivo infection is less understood. METHODS: Here, we developed an efficient SARS-CoV-2-RBD binding assay suitable for super resolution microscopy and simultaneous hACE2 immunodetection and mapped the correlation between hACE2 receptor abundance and SARS-CoV-2-RBD binding, both in vitro and in human lung biopsies. Next, we explored the specific proteome of SARS-CoV-2-RBD/hACE2 through a comparative mass spectrometry approach. FINDINGS: We found that only a minority of hACE2 positive spots are actually SARS-CoV-2-RBD binding sites, and that the relationship between SARS-CoV-2-RBD binding and hACE2 presence is variable, suggesting the existence of additional factors. Indeed, we found several interactors that are involved in receptor localization and viral entry and characterized one of them: SLC1A5, an amino acid transporter. High-resolution receptor-binding studies showed that co-expression of membrane-bound SLC1A5 with hACE2 predicted SARS-CoV-2 binding and entry better than hACE2 expression alone. SLC1A5 depletion reduces SARS-CoV-2 binding and entry. Notably, the Omicron variant is more efficient in binding hACE2 sites, but equally sensitive to SLC1A5 downregulation. INTERPRETATION: We propose a method for mapping functional SARS-CoV-2 receptors in vivo. We confirm the existence of hACE2 co-factors that may contribute to differential sensitivity of cells to infection. FUNDING: This work was supported by an unrestricted grant from "Fondazione Romeo ed Enrica Invernizzi" to Stefano Biffo and by AIRC under MFAG 2021 - ID. 26178 project - P.I. Manfrini Nicola.

Host range and structural analysis of bat-origin RshSTT182/200 coronavirus binding to human ACE2 and its animal orthologs.

Bat-origin RshSTT182 and RshSTT200 coronaviruses (CoV) from Rhinolophus shameli in Southeast Asia (Cambodia) share 92.6% whole-genome identity with SARS-CoV-2 and show identical receptor-binding domains (RBDs). In this study, we determined the structure of the RshSTT182/200 receptor binding domain (RBD) in complex with human angiotensin-converting enzyme 2 (hACE2) and identified the key residues that influence receptor binding. The binding of the RshSTT182/200 RBD to ACE2 orthologs from 39 animal species, including 18 bat species, was used to evaluate its host range. The RshSTT182/200 RBD broadly recognized 21 of 39 ACE2 orthologs, although its binding affinities for the orthologs were weaker than those of the RBD of SARS-CoV-2. Furthermore, RshSTT182 pseudovirus could utilize human, fox, and Rhinolophus affinis ACE2 receptors for cell entry. Moreover, we found that SARS-CoV-2 induces cross-neutralizing antibodies against RshSTT182 pseudovirus. Taken together, these findings indicate that RshSTT182/200 can potentially infect susceptible animals, but requires further evolution to obtain strong interspecies transmission abilities like SARS-CoV-2.

Mass photometry reveals SARS-CoV-2 spike stabilisation to impede ACE2 binding through altered conformational dynamics.

Here we show using mass photometry how proline substitutions, commonly used for SARS-CoV-2 spike stabilisation in vaccine design, directly affects ACE2 receptor interactions via dynamics of open and closed states. Conformational changes and ACE2 binding were influenced by spike variant and temperature, but independent of site-specific N-glycosylation.

HIV-1 gp120-CXCR4 recognition probed with synthetic nanomolar affinity D-peptides containing fragments of gp120 V3 loop.

The human immunodeficiency virus type 1 (HIV-1) recognizes one of its principal coreceptors, the CXC chemokine receptor 4 (CXCR4) on the host cell via the third variable loop (V3 loop) of HIV-1 envelope glycoprotein gp120 during the viral entry process. Here, we investigated the stereochemical mechanism of the molecular recognition of HIV-1 gp120 V3 loop with coreceptor CXCR4 by using peptide probes containing important fragments of the V3 loop. The tip and base/stem fragments of the V3 loop critical for V3 loop function were linked individually with the fragment derived from another CXCR4's chemokine ligand, vMIP-II to generate nanomolar affinity peptide probes of the interactions of CXCR4-V3 loop fragments. When the amino acid residues of the V3 loop fragments in these combinational peptides were changed from L-to D-configurations, the resulting peptides remarkably retained or had even enhanced recognition by CXCR4 as shown by competitive ligand-receptor binding. The ability of these peptides, regardless of the different l- or d-amino acids used, in binding CXCR4 and antagonizing CXCR4 functions was demonstrated by their blockade of calcium influx, cell migration, and CXCR4 internalization triggered by the activation of CXCR4 signaling by its endogenous ligand SDF-1α. The structural mechanisms of CXCR4 interactions with these peptides were examined with site-directed mutagenesis and molecular modeling. These results indicate that CXCR4's interface with key segments of HIV-1 gp120 V3 loop is flexible in terms of stereospecificity of ligand-receptor interaction which may have implication on understanding the viral entry mechanism and how the virus evades immune detection with V3 loop mutations and retains effective recognition of the host cell's coreceptor.

Biophysical and structural characterizations of the effects of mutations on the structure-activity relationships of SARS-CoV-2 spike protein.

Mutations on the spike (S) protein of SARS-CoV-2 could induce structural changes that help increase viral transmissibility and enhance resistance to antibody neutralization. Here, we report a robust workflow to prepare recombinant S protein variants and its host receptor angiotensin-convert enzyme 2 (ACE2) by using a mammalian cell expression system. The functional states of the S protein variants are investigated by cryo-electron microscopy (cryo-EM) and negative staining electron microscopy (NSEM) to visualize their molecular structures in response to mutations, receptor binding, antibody binding, and environmental changes. The folding stabilities of the S protein variants can be deduced from morphological changes based on NSEM imaging analysis. Differential scanning calorimetry provides thermodynamic information to complement NSEM. Impacts of the mutations on host receptor binding and antibody neutralization are in vitro by kinetic binding analyses in addition to atomic insights gleaned from cryo-electron microscopy (cryo-EM). This experimental strategy is generally applicable to studying the molecular basis of host-pathogen interactions.

Omicron BA.1 and BA.2 variants increase the interactions of SARS-CoV-2 spike glycoprotein with ACE2.

SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. First the BA.1 and later the BA.2 variant became the most dominant strains of the Omicron variant. To investigate how the mutations of these strains affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water, and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ∼34% and ∼51% stronger binding free energies to PD for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells.

Ligand-based design of peptide entry inhibitors targeting the endosomal receptor binding site of filoviruses.

Filoviruses enter cells through macropinocytosis and trafficking into the endosomes in which they bind to the receptor Niemann-Pick C1 protein (NPC1) for membrane fusion and entry into the cytoplasm. The endosomal receptor-binding is critical step for filovirus entry. Designing inhibitors to block receptor binding will prevent viral entry. Using available binding structural information from the co-crystal structures of the viral GP with the receptor NPC1 or with monoclonal antibodies, we have conducted structure-based design of peptide inhibitors to target the receptor binding site (RBS). The designed peptides were tested for their inhibition activity against pseudo-typed or replication-competent viruses in a cell-based assay. The results indicate that these peptides exhibited strong activities against both Ebola and Marburg virus infection. It is expected that these peptides can be further developed for therapeutic use to treat filovirus infection and combat the outbreaks.

Influenza A Virus Agnostic Receptor Tropism Revealed Using a Novel Biological System with Terminal Sialic Acid Knockout Cells.

Avian or human influenza A viruses bind preferentially to avian- or human-type sialic acid receptors, respectively, indicating that receptor tropism is an important factor for determining the viral host range. However, there are currently no reliable methods for analyzing receptor tropism biologically under physiological conditions. In this study, we established a novel system using MDCK cells with avian- or human-type sialic acid receptors and with both sialic acid receptors knocked out (KO). When we examined the replication of human and avian influenza viruses in these KO cells, we observed unique viral receptor tropism that could not be detected using a conventional solid-phase sialylglycan binding assay, which directly assesses physical binding between the virus and sialic acids. Furthermore, we serially passaged an engineered avian-derived H4N5 influenza virus, whose PB2 gene was deleted, in avian-type receptor KO cells stably expressing PB2 to select a mutant with enhanced replication in KO cells; however, its binding to human-type sialylglycan was undetectable using the solid-phase binding assay. These data indicate that a panel of sialic acid receptor KO cells could be a useful tool for determining the biological receptor tropism of influenza A viruses. Moreover, the PB2KO virus experimental system could help to safely and efficiently identify the mutations required for avian influenza viruses to adapt to human cells that could trigger a new influenza pandemic. IMPORTANCE The acquisition of mutations that allow avian influenza A virus hemagglutinins to recognize human-type receptors is mandatory for the transmission of avian viruses to humans, which could lead to a pandemic. In this study, we established a novel system using a set of genetically engineered MDCK cells with knocked out sialic acid receptors to biologically evaluate the receptor tropism for influenza A viruses. Using this system, we observed unique receptor tropism in several virus strains that was undetectable using conventional solid-phase binding assays that measure physical binding between the virus and artificially synthesized sialylglycans. This study contributes to elucidation of the relationship between the physical binding of virus and receptor and viral infectivity. Furthermore, the system using sialic acid knockout cells could provide a useful tool to explore the sialic acid-independent entry mechanism. In addition, our system could be safely used to identify mutations that could acquire human-type receptor tropism.

Characterizing Receptor Flexibility to Predict Mutations That Lead to Human Adaptation of Influenza Hemagglutinin.

A key step in the emergence of human pandemic influenza strains has been a switch in binding preference of the viral glycoprotein hemagglutinin (HA) from avian to human sialic acid (SA) receptors. The conformation of the bound SA varies substantially with HA sequence, and crystallographic evidence suggests that the bound SA is flexible, making it difficult to predict which mutations are responsible for changing HA-binding preference. We performed molecular dynamics (MD) simulations of SA analogues binding to various HAs and observed a dynamic equilibrium among structurally diverse receptor conformations, including conformations that have not been experimentally observed. Using one such novel conformation, we predicted─and experimentally confirmed─a set of mutations that substantially increased an HA's affinity for a human SA analogue. This prediction could not have been inferred from the existing crystal structures, suggesting that MD-generated HA-SA conformational ensembles could help researchers predict human-adaptive mutations, aiding surveillance of emerging pandemic threats.